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Abstract 


Vibrational  excitation  and  relaxation  for  a  molecule  adsorbed  on  a  surface 
is  investigated  in  connection  with  desorption/dissociation  processes.  The 
population  equations  for  a  laser-driven  anharmonic  oscillator  involving  both 
energy-  and  phase-relaxation  mechanisms  are  derived  with  the  Zwanzig  projector 
technique  and  eigenfunction-expansion  method  due  to  Weidlich.  Phase  relaxation 
is  seen  to  assist  in  the  excitation  of  the  resonant  active  mode. 


I .  Introduction 

Infrared  laser  radiation  impinging  onto  a  molecule  adsorbed  on  a  solid 
surface  can  resonantly  excite  the  admolecular  internal  vibrational /rotational 
modes,  or  deposit  photon  energy  in  the  adbond  with  which  the  molecule  is  bound 
to  the  surface. It  has  been  shown  experimentally^’^  that  SFg  and  CH^F 
molecules  adsorbed  on  NaCl  surfaces  at  low  temperatures  can  be  desorbed  by 
the  resonant  excitation  of  the  adsorbate  internal  vibration,  namely,  the 
mode.  The  treatment  of  the  problem  of  many  degrees  of  freedom  for  the 
molecule-surface  system  is  very  complicated  and  generally  an  insoluble  problem. 
Interactions  between  the  resonant  active  vibrational  mode  and  the  other  inactive 
modes,  and  between  the  admolecule  and  surface  excitations  -  surface  phonons,  surface 
plasmons,  reflected  photons  (i.e.,  photons  emitted  by  the  admolecule  and 
reflected  by  the  surface),  etc.  -  create  a  large  number  of  dissipative  channels 
by  which  absorption,  dissociation,  desorption  or  chemical  reaction  on  the  surface 
are  strongly  influenced.  (Desorption  is  the  breaking  of  the  adbond  between  the 
admolecule  and  the  surface,  and  dissociation  is  the  breaking  of  a  bond  within 
the  admolecule.)  In  this  Letter,  we  are  especially  interested  in  the  influence 
of  surface-dephasing  relaxation  on  the  excitation  of  the  resonant  mode,  which 
plays  an  important  role  in  the  process  of  desorption/dissociation. 


Theory 


Under  normal  circumstances,  due  to  the  vast  number  of  degrees  of  freedom 
of  the  admolecul e-surface  system,  the  intrasystem  relaxation  can  be  very  fast,  and  in 
fact,  comparable  to  or  faster  than  the  rate  of  photon  absorption.  In  this  case, 
energy  randomization  or  statistical  thermodynamical  behavior  may  at  least  partly 

Q 

occur  in  desorption/dissociation  processes  on  surfaces.  The  IR-laser-driven 
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resonant  active  mode  absorbs  photons  and  transfers  its  energy  to  other 
vibrational  modes  or  the  surface,  and  then  "interference"  superposition  of 
these  many  vibrations  can  instantaneously  accumulate  energy  in  the  weakest  bond 
to  lead  to  molecule  dissociation  or  desorption  from  the  surface,  if  this  energy 


is  beyond  the  appropriate  theshold 


9,10 


As  an  example,  the  threshold  temperature  for 

.  i  i  _  i  «  ly  B 


thermal  dissociation  of  SFg  in  the  gas  phase  is  measured  to  be  1600  K. 

At  this  temperature,  the  average  occupation  number  of  the  resonant  mode  (v-j) 

of  SFC  is  n=0.76,  which  means  that,  while  the  energy  deposited  in  the  molecule 
0 

has  gone  into  the  mode  to  exceed  the  critical  value  for  dissociation,  the 


other  modes  remain  in  low  excited  states. 


We  shall  confine  ourselves  here  to  the  excitation  and  relaxation  of  the 
active  anharmonic  vibrational  mode  which  is  driven  by  a  resonant  coherent  field 
and  controlled  by  both  energy  (Tj)  and  phase  (T^)  relaxation.  For  this  purpose, 
we  treat  the  many  other  non-active  modes  and  surface  excitations  (i.e.,  phonons, 
plasmons,  reflected  photons,  etc.)  as  a  large  thermal  reservoir,  by  which 
dissipation  mechanisms  are  provided.  A  reduced  density  operator  W  for  the 
resonant  mode  obeys  the  master  equation 


dW  s  fdw"j  +  fdw"[ 

dt  Ht  coherent  MM1 


ncoherent 


■ilc  +  Aj  +  Ag  W. 


The  reversible  term  is  -iLcW  =  -(i/tf)[H,W],  where  H  =  HM  +  Hp  is  the  sum  of 


anharmonic  oscillator  and  external  field  interactions. 


Hr  *  )a  a  -  'tTea~a(a  a+1) 


Hp  =  wE(t)(a  +a). 


(3) 


Here  a*  and  a  are  oscillator  ladder  operators,  y  is  the  molecular  dipole  moment, 

E(t)  is  the  applied  electric  field,  is  the  laser  frequency,  and  e  is  a 

12  13 

positive  anharmonicity.  The  dissipative  terms  are  given  by  * 


AjW  *  y+([a,Wa^]  +  [aW,a  ]) 

+  YA{[a~,Wa]  +  [a  W,a]) ,  (4) 

A«W  *  n([a~a,Wa  a]  +  [a  aW,  a  a]).  (5) 

14 

AgW  is  responsible  for  a  pure  dephasing  in  the  decay  of  the  off-diagonal  matrix  elements, 
and  yf,  y+  and  n  contain  implicit  information  about  the  reservoir,  including  the  surface. 
With  the  dissipative  terms  [Eqs.  (4)  and  (5)  only],  the  relaxation  behavior  of  the 
expectation  values  of  the  amplitude  and  number  operator  is  governed  by  the  equations 
m<a>s-(Y,-Y  A+n)<a>  and  -^<a  a> = 2(Y>-Y*)<a+a> + 2yf.  A2  thus  affects  only  the  off-diagonal 
elements  of  the  density  operator  by  providing  a  pure  phase-destroying  term  with  a  rate 
of  n.  In  the  Born  approximation,  the  diagonal  parts  of  the  density  matrix  obey  the  equation 


where  is  a  Liouville  operator  due  to  the  oscillator-field  interaction,  P  is  a 

projection  operator  used  to  obtain  the  diagonal  part  of  the  density  operator,15 

and  Wr  =  PW  is  the  relevant  part  of  the  density  matrix. 

The  difficult  problem  here  is  the  explicit  evaluation  of  the  matrix  elements 

that  appear  inside  the  Integral  In  Eq.  (6).  We  use  the  Weld  1 ich  elgenfuntlon- 
12  13 

expansion  method  ’  to  evaluate  these  elements.  The  "eigenstate"  of  the  super¬ 
operators  a2  and  A2  are  I !APq)  and  i iBp^), respectively,  and  under  the  Markoff 
assumption,  the  "eigenvalues"  of  A2  have  the  form  xpg  -  -nq*,  where  p  ■  0,  1,  2,  . 
and  q  »0,  ±1,  ±2,  •••.  The  equations  of  motion  for  the  diagonal  parts  of  the 
density  matrix  are  then  given  as 


5 


Wn  =  <n).'.1Wr|n>  -  dt  l  e"^ 

J0  pq 

*  Tr(Bpq  LiWr(t))  <niLie  X^Apq ln>*  (V 

By  a  lengthy  but  straightforward  algebraic  calculation,  we  obtain  a  simpler  set 
of  coupled  differential  equations  to  describe  the  excitation  and  relaxation 
behavior  for  the  resonant  vibrational  mode: 


3p(t)  , 

2*  C(n+1)  W„+1  -  nHn]  -  2r,  [-=—  ]2 


<n*‘> 


f  "(VVl> 
t1  *  t^<n-l>  *  § 


1  +  (2  ^  n  -  ±)‘ 

n  n 


The  decay  parameters  <  and  n  are  determined  by  the  internal  mode-mode 

coupling  and  the  admolecule-surface  interaction,  and  for  a  fixed  temperature  and 
admolecule-surface  distance  they  can  be  considered  as  constants  (here  a  cold 
surface  of  0  K  has  been  assumed),  n^(t)  is  the  Rabi  frequency,  and  A  is  the 
detuning  between  the  admolecular  fundamental  and  laser  frequencies. 

The  first  two  terms  with  the  factor  of  2<  on  the  right-hand  side  of  Eq.  (8) 
describe  the  energy  relaxation  processes  which  always  lead  to  population  loss 
in  the  resonant  mode.  In  the  next  terms  with  the  factors  of  n  describing  the 
influence  of  the  phase  relaxation,  we  notice  that  there  are  competitive  factors: 
upon  an  Increase  of  the  value  of  n.  the  phase  relaxation  tends  to  diminish  the 
detuning  and  anharmonicity  "bottleneck”  effects,  and  the  effective  Rabi  frequency 
f;R/n*  tends  to  decrease. 


Ill .  Calculations  and  Results 

Figure  1  shows  the  influences  of  the  phase  relaxation.  The  time  evolution 
for  the  average  occupation  number  n  of  a  driven  damped  anharmonic  oscillator  is 
displayed,  with  consecutive  values  of  the  phase-relaxation  parameter  n  (n= 

1,  5,  9,  13).  We  note  that  there  is  a  range  of  values  of  n  which  assist  the 
excitation  of  the  admolecular  active  mode.  For  the  optimal  parameter  ns5  in 
our  specific  case,  Figure  2  displays  the  time-dependent  population  distributions 
for  several  excited  states,  where  considerable  population  is  seen.  One  would 
generally  expect  the  phase  decay  r\  to  be  much  more  enhanced  than  k  by  the 
inactive  adspecies  vibrational  modes  and  by  the  surface.^  The  phase-dissipation 
mechanism  may  be  generated  in  several  ways  in. an  admolecul e-surface  system  as 
follows:  (a)  intramolecular  coupling  between  the  active  and  inactive  vibrational 
modes; 17  (b)  intermolecular  collisions  when  the  admolecules  can  migrate  on  the 
surface;  (c)  rotational  relaxation  within  the  vibrational  level  manifold  of 

the  admolecule;^  (d)  self-coupling  of  the  admolecule  due  to  surface- reflected 
photons.  19  >20 
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1.  Time  dependence  of  the  average  occupation  number  rf  =  <a+a>  for  (cw)  laser- 
driven  anharmonic  oscillator  in  contact  with  a  bath.  The  curves  correspond 
to  n  =  1,  5,  9,  13,  and  the  others  parameters  are  (<,  fiR,  e.  A)  =  (0.01, 

13 

0.5,  1.5,  3.0),  where  all  rates  are  in  the  unit  of  the  transverse  relaxation 
rate.  This  figure  shows  that  when  the  vibrational  linewidth  n  is  smaller 
than  the  detuning  A,  the  maximum  vibrational  excitation  is  small  but  grows 
as  n  increases.  Excessively  large  values  of  n  eventually  reduce  the  average 
excitation  level,  presumably  due  to  the  spread  of  population  within  the  line- 
width. 

2.  The  probability  of  occupation  of  level  n,  Wp ( t ) ,  as  a  function  of  time,  n 
is  fixed  at  5,  and  the  other  parameters  are  given  in  Figure  1. 
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